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ABSTRACT 

This paper presents a comprehensive examination of the X-ray continuum variability of the 
bright Seyfert 1 galaxy MCG-6-30-15. The source clearly shows the strong, linear correlation 
between rms variability amplitude and flux first seen in Galactic X-ray binaries. The high 
frequency power spectral density (PSD) of MCG-6-30-15 is examined in detail using a Monte 
Carlo fitting procedure and is found to be well represented by a steep power-law at high 
frequencies (with a power-law index a « 2.5), breaking to a flatter slope (a ss 1) below 
/br ~ 0.6 — 2.0 X 10“"^ Hz, consistent with the previous results of Uttley, M'^Hardy & 
Papadakis. The slope of the power spectrum above the break is energy dependent, with the 
higher energies showing a flatter PSD. At low frequencies the variations between different 
energy bands are highly coherent while at high frequencies the coherence is significantly 
reduced. Time lags are detected between energy bands, with the soft variations leading the 
hard. The magnitude of the lag is small (^ 200 s for the frequencies observed) and is most 
likely frequency dependent. These properties are remarkably similar to the temporal properties 
of the Galactic black hole candidate Cygnus X-1. The characteristic timescales in these two 
types of source differ by ~ 10^; assuming that these timescales scale linearly with black 
hole mass then suggests a black hole mass ^ 10® Mq for MCG-6-30-15. We speculate that 
the timing properties of MCG-6-30-15 may be analogous to those of Cyg X-1 in its high/soft 
state and discuss a simple phenomenological model, originally developed to explain the timing 
properties of Cyg X-1, that can explain many of the observed properties of MCG-6-30-15. 
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15 - X-ray: galaxies 


1 INTRODUCTION 

The X-ray emission from radio-quiet Active Galactic Nuclei 
(AGN) shows persistent variability that is stronger and more rapid 
than in any other waveband. This rapid X-ray variability provided 
early support for the black hole/accretion disc model of AGN (e.g. 
Rees 1984) and suggested that the X-ray emission is produced in 
the inner regions of the central engine (see Mushotzky, Done & 
Pounds 1993 for a review). The X-ray fluctuations appear erratic 
and random, with variability occurring on a range of timescales 
(e.g. M'^Hardy 1989). 

In the study of time-variable processes one of the most widely 
used statistical tools is the variability power spectrum (e.g. Priest¬ 
ley 1981; Bloomfield 2000). The power spectral density (PSD) de¬ 
scribes the amount of variability power (i.e. amplitude^) as a func¬ 
tion of temporal frequency (timescale”^). The long, uninterrupted 
observations provided by EXOSAT allowed for the measurement 
of the first accurate X-ray PSDs for Seyfert 1 galaxies. These EX¬ 
OSAT PSDs were shown to rise at lower frequencies as a power- 
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law: 'P{f) oc /”“, where Vif) is the power at frequency / and a 
is the PSD slope, found to be a ~ 1.5 over the frequencies ~ 10“® 
to ~ 10“® Hz (Green, M'^Hardy & Lehto 1993; Lawrence & Pa¬ 
padakis 1993). This type of broad-band variability, rising to lower 
frequencies with a power spectral slope a ^ 1, is usually called 
“red noise” (for an introduction to red noise in astronomy and else¬ 
where see Press 1978). 

The red noise PSDs of Seyfert Is are similar to those observed 
in Galactic Black Hole Candidates (GBHCs, see e.g. van der Klis 
1995) on much shorter timescales, and suggests the X-ray emis¬ 
sion mechanisms may be the same in these objects that differ in 
black hole mass by factors of ^ 10® (M'^Hardy 1989). The PSD of 
the best studied GBHC, Cygnus X-1 in its low/hard state, can be 
approximated by a doubly-broken power-law with a steep slope at 
high frequencies (a « 2 above ~ 3 Hz), breaking to a flatter slope 
at intermediate frequencies (a « 1 in the range ~ 0.2 — 3 Hz) then 
breaking to a flat (“white”) spectrum at low frequencies (a « 0 
below 0.2 Hz). See Nowak et al. (1999a) and Belloni & Hasinger 
(1990) for more details of Cyg X-1. Other GBHCs such as GX 339- 
4 show similar behaviour (Nowak, Wilms & Dove 2002). The posi¬ 
tion of the breaks in the PSD represent “characteristic timescales” 
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of the system, but these timescales have no clear physical inter¬ 
pretation at present (e.g. Belloni, Psaltis & van der Klis 2002 and 
references therein). 

The featureless iiXOSArPSDs of Seyfert Is provided no such 
timescales. Such breaks must exist, however, or the total power 
would diverge at low frequencies (for PSD slopes a ^ 1). In¬ 
deed, in recent years dedicated X-ray monitoring programmes with 
RXTE have revealed deviations from the power-law PSD in a few 
Seyfert 1 galaxies; at low frequencies the PSD breaks to a flat¬ 
ter slope (Edelson & Nandra 1999; Uttley, M'^Hardy & Papadakis 
2002; Markowitz et al. 2002). The timescale of this break fre¬ 
quency is consistent with the idea that the PSD of AGN and GB- 
HCs are essentially the same, with the characteristic timescales 
scaling linearly with the black hole mass. Such a linear scaling of 
timescales with black hole mass might be expected if the character¬ 
istic timescales are related to e.g. a light-crossing timescale or an 
orbital timescale. This linear scaling also means that a 100 ks light 
curve of a Seyfert galaxy (such as presented here) is comparable to 
a < 1 s light curve of Cyg X-1. On these timescales a bright Seyfert 
1 galaxy provides many more photons (by more than two orders of 
magnitude) than the GBHC, i.e. Seyferts have a higher count rate 
per characteristic timescale. This means that intensive monitoring 
of Seyferts can provide a view of the very high frequency variations 
in accreting black holes that is difficult to access even in bright GB- 
HCs. 

Not only do AGN and GBHCs show similar PSDs, but they 
display other striking similarities in their X-ray timing proper¬ 
ties. Uttley & M'^Hardy (2001) showed that both X-ray binary and 
AGN light curves show a tight, linear correlation between flux and 
rms variability amplitude (see also Edelson et al. 2002). The vari¬ 
ability properties of GBHCs such as Cyg X-1 are a function of 
photon energy. The PSD shape above the high frequency break 
is energy-dependent, with the harder energy bands showing flat¬ 
ter PSD slopes (e.g. Nowak et al. 1999a; Lin et al. 2000). In ad¬ 
dition the variations between different bands are highly coherent 
(well correlated at each Fourier frequency) and show time delays, 
with variations at softer energies leading those at harder energies 
(Miyamoto & Kitamoto 1989; Cui et al. 1997; Nowak et al. 1999a). 
The magnitude of the time delay is frequency dependent, with the 
delay increasing at lower frequencies. To date these properties have 
only been measured in one AGN, NGC 7469 (Nandra & Papadakis 
2001; Papadakis, Nandra & Kazanas 2001), and in this object the 
energy-dependence of the PSD slope, the coherence and time de¬ 
lay properties appeared remarkably similar to those seen in Cyg 
X-1. These similarities re-enforce the idea that the X-ray emission 
mechanisms are the same in GBHCs and AGN. 

This paper presents a detailed analysis of the X-ray contin¬ 
uum variability of the bright, variable Seyfert 1 galaxy MCG-6- 
30-15 (z = 0.007749) using a long XMM-Newton observation. 
The large collecting area, wide band-pass and long orbit (~ 2 day) 
of XMM-Newton make it ideal for examining the high frequency 
(^ 10~® Hz) timing properties of Seyfert galaxies. The analysis 
concentrates on frequency domain (i.e. Fourier) time series analysis 
techniques, as these are used extensively in the analysis of GBHC 
data. The rest of this paper is organised as follows. In Section ^ 
the data reduction procedures are outlined. In Section ^ the basic 
properties of the light curves are discussed, followed by an analy¬ 
sis of the PSD in Section Spectral variability is then examined 
in detail by measuring the coherence and time lags between dif¬ 
ferent energy bands, as discussed in Section As these analyses 
are relatively new to AGN studies we discuss the methodology in 


some detail. Finally, the implications of these results are discussed 
in Section 1^ and some concluding remarks are given in Section^ 


2 DATA REDUCTION 

XMM-Newton (Jansen et al. 2001) observed MCG-6-30-15 over 
the period 2001 July 31 - 2001 August 5 (revolutions 301, 302 
and 303), during which all instruments were operating nominally. 
Fabian et al. (2002) discuss details of the observation and basic 
data reduction. For the present paper the data were reprocessed en¬ 
tirely with the latest software (SAS v5.3.3) but the procedure 
followed that discussed in the previous paper (except that the ex¬ 
traction region was a circle of radius 35 arcsec for both MOS and 
pn). 

Light curves were extracted from the EPIC pn data in four dif¬ 
ferent energy bands: 0.2-10.0 keV (full band), 0.2-0.7 keV (soft 
band), 0.7-2.0 keV (medium band) and 2.0-10.0 keV (hard band). 
These were corrected for telemetry drop outs (less than 1 per cent 
of the total time), background subtracted and binned to 100 s time 
resolution. The errors on the light curves were calculated by propa¬ 
gating the Poisson noise. The light curves were not corrected for the 
~ 71 per cent “live time” of the pn camera (Struder et al. 2001), 
which is only a scaling factor. The EPIC MOS light curves were 
also examined, but these are essentially identical to the pn light 
curves except with lower signal-to-noise. Therefore, in the rest of 
this paper the analysis concentrates on the pn data, but the MOS 
data are used to confirm the pn results of Section The full band 
light curves are shown in Fig. The revolutions provided three 
uninterrupted light curves of 80, 123 and 124 ks duration. 

As discussed in Fabian et al. (2002), BeppoSAX observed 
MCG-6-30-15 simultaneously with XMM-Newton, but these data 
are not used in the present paper. The LEGS and MECS data cover 
approximately the same bandpass as the EPIC (but with lower sig¬ 
nal to noise) and the higher energy PDS data (with a count rate 
~ 0.7 ct s“^) do not have sufficient statistics to determine the high 
energy variability properties accurately. In addition, analysis of the 
BeppoSAX data is complicated by the poor sampling of the light 
curves (due to repeated Earth occulatations and SAA passages), in¬ 
deed the amount of “good” exposure time obtained for the PDS was 
40 ksec spread over a > 400 ksec observing period. 


3 BASIC TEMPORAL PROPERTIES 

The source shows a factor of five variation between minimum and 
maximum flux (both occurring during revolution 303), and signifi¬ 
cant variability can be seen on timescales as short as ~ 100 s. 

As a first measure of the variability amplitude the root mean 
square (rms) amplitude, in excess of Poisson noise, was calculated 
as 

O’rms — \/ rTerr, ( 1 ) 

where is the total variance of the light curve x(ti) of length N 
data points 

N 

^ - X)^, ( 2 ) 

i=l 

where x is the mean of x{ti) and alri. is the contribution expected 
from measurement errors 
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Figure 1. Full-band (0.2-10.0 keV) EPIC pn light curves in 100 s bins. The top panel shoves all three consecutive light curves and the bottom panels sho3v the 
individual light curves from each XMM-Newton revolution. 


N 

O-err = jy X! (^i)l^ ’ (3^ 

i=l 

and CTerr(ii) IS the error on x{ti). 

The fractional excess rms variability amplitude (TVar = 
Crms/s; Edelson et al. 2002) and mean count rate (x) were calcu¬ 
lated for each of the three revolutions and are tabulated in Table [i[ 
The errors on TVar were estimated using the prescription given by 
Edelson et al. (2002; but see their appendix A for caveats on the 
interpretation of these estimates). The medium band showed the 
highest variability amplitude. This can also been seen in Fig. 6 of 
Fabian et al. (2002), where the fractional rms spectrum peaked in 
the 0.7-2.0 keV band. 


3.1 Radiative efficiency 

Fig. I shows a rapid “event” from the revolution 302 light curve. 
During this the source changed in 0.2-10 keV luminosity by AL « 
1.5 X 10“^^ erg s“^ in 100 s. This was used to estimate the radiative 
efficiency rj of the source, assuming photon diffusion through a 
spherical mass of accreting matter: AL/At < rj ■ 2.1 x 10“*^ erg 
s“^ (Fabian 1979; Brandt et al. 1999). This gave a limit of 77 ^ 0.7 
per cent, below the expected maximum for accretion onto a black 
hole (6 and 30 per cent for non-rotating and maximally rotating 
black holes, respectively; Thome 1974). 


Table 1. Mean count rates (x) and fractional excess rms variability ampli¬ 
tudes (Fvar) for the light curves from each revolution and in each energy 
band. 


Revolution 

band 

X (ct s 

Uvar (%) 

301 

Full 

18.51 

28.7 ±0.7 


Soft 

8.63 

27.8 ± 0.7 


Medium 

6.68 

32.3 ± 0.8 


Hard 

3.27 

26.0 ±0.7 

302 

Full 

19.94 

19.3 ± 0.4 


Soft 

8.85 

18.4 ±0.4 


Medium 

7.53 

21.4 ±0.4 


Hard 

3.61 

19.3 ± 0.4 

303 

Full 

18.41 

34.9 ±0.7 


Soft 

8.56 

31.8 ±0.6 


Medium 

6.69 

40.1 ± 0.8 


Hard 

3.22 

33.7 ±0.7 


3.2 rms-fiux correlation 

Using the 100 s resolution light curves the mean count rate and 
(Trms were calculated in bins of 15 data points. The values of arms 
from all three orbits were then binned as a function of count rate 
such that there were 20 measurements per bin, and an error was as- 
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Figure 2. Section of the full-band (0.2-10.0 keV) EPIC pn light curve 
(from rev. 302) showing a rapid “flare” event. 
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Figure 3. Correlation between nns variability amplitude and flux. 


Table 2. Results of fitting a linear function to the rms-flux correlation. Er¬ 
rors on the model pai'ameters con'espond to a 90 per cent confidence level 
for one interesting parameter (i.e. a Ax^ = 2.7 criterion). 


Band 

y-intercept 

Gradient 

X^/dof 

Full 

0.07 ±0.22 

0.052 ±0.013 

10.66/10 

Soft 

0.02 ±0.12 

0.049 ±0.015 

7.78/10 

Medium 

0.04 ±0.07 

0.057 ±0.016 

6.97/10 

Hard 

0.001 ± 0.049 

0.065 ±0.016 

6.68/10 


signed based on the scatter within each bin (equation 4.14 of Bev- 
ington & Robinson 1992). Fig. ^ shows the resulting correlation 
between arms and count rate in the four energy bands. The signifi¬ 
cance of the correlation was tested using the Spearman rank-order 
correlation coefficient, and the Kendall r coefficient (Press et 
al. 1992). Both these indicate that the correlation is significant at 
^ 99.5 per cent confidence in all energy bands. The effect of the 
rms-flux correlation can be seen in the light curves (Fig. 0): the 
peaks of the light curves appear more variable (“jagged”) while the 
troughs appear much smoother. 

Table ^ shows the results of fitting a linear model to these 
data; in all cases the linear model gives a good fit to the data, with 
a y-axis intercept consistent with zero. The gradient of this relation 
corresponds to the TVar measured only on these timescales (be¬ 
tween 100 and 1500 s). The la upper limit on any constant flux 
offset (i-intercept) is 1.2 ct s“^ for the full band, or 6.2 per cent of 
the total count rate in that band. In the three energy sub-bands the 
upper limits on the constant flux offset are as follows: 1.04 ct s“^ 
(12 per cent), 0.09 ct s“^ (1.3 per cent) and 0.43 ct s“^ (12.8 per 
cent) for the soft, medium and hard bands respectively. 


4 POWER SPECTRAL PROPERTIES 
4.1 Estimating the power spectrum 

The PSD (sometimes known as the auto spectrum) describes the 
amount of variability “power” as a function of Fourier frequency. 
The PSD and auto correlation function (ACF) are Fourier pairs, i.e. 
they are Fourier transforms of one another. For a general introduc¬ 
tion to Fourier methods of time series analysis see e.g. Priestley 


(1981) and Bloomfield (2000), and for a review of Fourier tech¬ 
niques as applied to X-ray time series analysis see van der Klis 
(1989). 

4.J.J Measuring the periodogram 

If the light curve x{ti) is evenly sampled (with a sampling pe¬ 
riod AT) then the PSD can be estimated by calculating the peri- 
odograrrfl. This is the normalised, modulus-squared of the Discrete 
Fourier Transform (DFT) of the data. The DFT of the light curve, 
X(fj), is calculated as (see Press et al. 1992): 

1 /o N 

( A rp K \_/2, 

— j = Rx{fj)+Ux{f,),W 

i=l 

where Rx{fj) and Ix{fj) are the real and imaginary parts of the 
DFT given by the discrete cosine and sine transforms, respectively: 

Rxifj) = (-^) cos{2nfjti) 

i=l 

= (—j '^x{ti) (5) 

i=l 

The transforms are normalised by (ATso that the am¬ 
plitudes are independent of the sampling and length of the light 
curve. The DFT is calculated at N/2 evenly spaced frequencies 
fi ~ j/NAT (where j = 1,2,..., N/2). Thus the highest fre¬ 
quency accessible to the DFT is the Nyquist frequency, /Nyq = 
/jv /2 = 1/2AT, and the lowest frequency is fi = 1/NAT. Note 
that it is customary to subtract the mean flux from the light curve 
before calculating the DFT, this eliminates the zero-frequency 
power. 

The complex valued DFT is then squared 
\X{fj)f = V*(/,)X(/,) = {Rx{f,)f + {Ix{m\ (6) 

^ Following Priestley (1981) we use the term “periodogram” to denote the 
discrete function P{fj) described by equation W which is merely an esti¬ 
mator of the continuous PSD Vif). The periodogram is therefore related 
to each new realisation of the process, whereas the PSD is representative of 
the true, underlying process. 
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where * denotes the complex conjugate. The periodogram, P{fj), 
is then calculated hy choosing an appropriate normalisation: 

P{ff = ^\Xifj)\\ (7) 

(Note that the Fourier transform has already been normalised in 
equations ^ and ^ to account for the data sampling.) 

This normalisation (defined in van der Klis 1997, see also 
Miyamoto et al. 1991) is most commonly used in analysis of AGN 
and X-ray binaries because the integrated periodogram yields the 
fractional variance of the data (i.e. This normalisation will 

be used throughout this paper. The factor of two is because this is a 
“one sided” normalisation; the periodogram is defined for both neg¬ 
ative and positive frequencies, with this normalisation the integral 
over only positive frequencies yields the fractional variance. The 
integrated periodogram between two frequencies, say /i and /2, 
yields the contribution to the fractional variance due to variations 
between the corresponding timescales I//2 and I//1 (this follows 
from ParsevaTs theorem, see e.g. van der Klis 1989). The units for 
the periodogram ordinate are (rms/mean)^ FIz“^ (where rms/mean 
is the dimensionless quantity TVar), or simply FIz“^. 


4.1.2 Poisson noise 

If the time series is a continuous photon counting signal binned into 
intervals of AT, as it is here, the effect of Poisson noise is to add 
an approximately constant amount of power to the periodogram 
at all frequencies. Thus, at high frequencies, where the red noise 
spectrum of the AGN is weak, the observed periodogram will be 
dominated by the flat (“white”) Poisson noise spectrum. With the 
above normalisation this constant Poisson noise level is Pn = 2/®. 
This formally assumes zero background flux, but for the present 
observation this is a reasonable approximation (even in the hard 
band the background count rate is ~ 0.6 per cent that of the source). 


4.1.3 Scatter in the periodogram 

A periodogram measured from a single light curve shows a great 
deal of scatter around the underlying PSD. This is a natural out¬ 
come of measuring the periodogram from a single realisation of 
a stochastic process. The result of this is that the periodogram at 
a given frequency, P{f), is scattered around the underlying PSD, 
'P(/), following a distribution with two degrees of freedom (see 
Section 6.2 of Priestley 1981): 

P{f) = P{f)xl/2, (8) 

where X 2 is a random variable distributed as x^ with two degrees of 
freedom, i.e. an exponential distribution with a mean and variance 
of two and four, respectivel5Q The expectation value of the peri¬ 
odogram at a frequency / is thus equal to the PSD but so is its stan¬ 
dard deviation, as a result the periodogram shows a great deal of 
scatter. This fundamental property of the periodogram is discussed 
terms of X-ray time series analysis by e.g. Leahy et al. (1983), van 
der Klis (1989), Papadakis & Lawrence (1993), Timmer & Konig 
(1995) and Stella et al. (1997). In addition the periodogram is an in¬ 
consistent estimator of the PSD because this intrinsic scatter does 


^ For the case of even N the Fourier transform at the Nyquist frequency 
(DTr(/Nyq)) is always real, and the periodogram at this frequency 
(P(/Nyq)) is distributed as xf, i-C- with one degree of freedom 


not decrease as the number of data points in the light curve is in¬ 
creased. The periodogram can, however, be averaged (binned over 
separate light curves or adjacent frequencies) to reduce this scat¬ 
ter and produce a consistent estimate of the PSD (Jenkins & Watts 
1968). 

In GBHC analysis the light curves are typically broken into 
many separate segments and periodograms are calculated for each 
segment. At a given Fourier frequency the periodogram estimates 
from each light curve are identically and independently distributed 
(assuming the variability is stationary) and so their average will be¬ 
come Gaussian (following the central limit theorem). Therefore, 
averaging the periodograms from many independent light curve 
segments will produce a consistent estimate of the PSD with Gaus¬ 
sian errors (e.g. van der Klis 1997). However, dividing up the light 
curve limits the lowest frequency probed by the periodogram and 
so is not so suitable for AGN analysis where it is crucial to examine 
the widest frequency range possible. 

An alternative technique is discussed by Papadakis & 
Lawrence (1993) and is in some sense “optimal” for red-noise 
light curves such as those from AGN. In this method a single pe¬ 
riodogram is computed from the light curve, reaching down to the 
lowest accessible frequencies, and periodogram estimates at con¬ 
secutive frequencies are averaged. For evenly sampled data the pe¬ 
riodogram estimates at each Fourier frequency are independently 
distributed (according to equation ^ and so the averaged peri¬ 
odogram estimate in each frequency bin will tend to a Gaussian 
distribution. In the method of Papadakis & Lawrence (1993) it is 
the logarithm of the periodogram that is binned, such that each 
bin contains A ^ 20 periodogram estimates, and an error is as¬ 
signed based on the scatter within each bin. Using the logarithm of 
the periodogram means the distribution of the binned periodogram 
converges on Gaussian with fewer data points per bin, and there¬ 
fore allows for finer frequency binning. This binning method is 
used throughout this paper to produce binned periodogram esti¬ 
mates with errors. (See also van der Klis 1997 for more on binned 
periodogram estimates.) 


4.1.4 Bias in the periodogram 

A further point is that periodograms measured from finite data tend 
to be biased by windowing effects which further complicate their 
interpretation (van der Klis 1989; Papadakis & Lawrence 1993; Ut- 
tley et al. 2002). The observed light curve (/(f)) is the true, con¬ 
tinuous light curve of the source (((f)) multiplied by the sampling, 
“window,” function (tu(f)): 

/(f) = ((f) • w{t), (9) 

where the window function takes the value unity when the light 
curve is sampled and zero elsewhere. From the convolution theo¬ 
rem of Fourier transforms, the Fourier transform of the observed 
light curve (T(/)) is the Fourier transform of the true light curve 
(!/(/)) convolved with the Fourier transform of the window func¬ 
tion (W{f)y. 

F{f) = Lif)<lDW{f). (10) 

The periodogram (which is the square of the discrete Fourier trans¬ 
form of the observed light curve) is therefore the periodogram of 
the continuous light curve of the source multiplied by the peri¬ 
odogram of the window function. Thus the periodogram distorted 
away from the true PSD by the sampling, i.e. the convolution with 
the spectral window function. 
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This leads to two main effects. The first is due to the discrete 
sampling of the data. If the light curve is not contiguously sampled 
then power from frequencies above the Nyquist frequency can be 
folded back, or “aliased,” into the observed frequency range. For 
binned, continuous light curves such as those presented here, the 
aliasing will negligible (see van der Klis 1989). In any event, from 
continuous data the very high frequency spectrum will dominated 
by the Poisson noise level and not aliasing. 

The second effect is “red noise leak” and arises as a result of 
the finite duration of the light curve. Power is transferred from low 
to high frequencies by the lobes of the spectral window function 
(see e.g. Deeter & Boynton 1982; van der Klis 1997). If there is 
significant power at frequencies below the lowest frequency probed 
by the periodogram (i.e. on timescales longer than the length of the 
observation), this can manifest itself as slow rising or falling trends 
across the light curve. These trends will contribute to the total vari¬ 
ance of the light curve (and thus the total power) but, as the peri¬ 
odogram does not extend to such low frequencies, the power will be 
transfered into the observed frequency band-pass (i.e. it will “leak” 
to higher frequencies). Red noise leak tends to add a constant com¬ 
ponent to the observed periodogram in the form of a power-law 
of slope a = 2 (i.e. the PSD of a linear trend; see van der Klis 
1997 for more details), the amplitude of this depends on the PSD 
of the source at low frequencies (see also section 3.3 of Papadakis 
& Lawrence 1995). The effect of red noise leak on the periodogram 
can be accounted for using the Monte Carlo simulation procedure 
outlined below. 


4.2 PSD fitting procedure 


The method employed in the present work is similar to the methods 
discussed by Done et al. (1992), Green, M'^Hardy & Done (1999) 
and Uttley et al. (2002). A model PSD is assumed, multiple realisa¬ 
tions (light curves) are simulated using this PSD and are resampled 
to match the original data (thereby including the distorting effects 
of the sampling). Periodograms are then calculated from the sim¬ 
ulated light curves in the same manner as the real data, averaged 
over the realisations, and this “average distorted model” (ADM), 
which includes the effects of red noise leak, is compared to the 
original periodogram using a goodness-of-fit measure. This method 
therefore accounts for the biases in the periodogram introduced by 
the window function. It is important to bear in mind that, as with 
any model-fitting procedure, the results depend on the assumptions 
made, i.e. the choice of specific models to test. 


The errors derived from the periodograms (see Section 4.1.3) 


are approximately Gaussian (Papadakis & Lawrence 1993) and 
therefore the statistic can be used to estimate the goodness of the 
fit. Once a reasonable fit is found, uncertainties on the model pa¬ 
rameters can be estimated using the Ax^ criterion (e.g. Lampton, 
Margon and Bowyer 1976), as in X-ray energy spectrum fitting. 

The algorithm described by Timmer & Kbnig (1995) was used 
to generate the simulated light curves from the model PSD. This 
algorithm generates a random time series from an arbitrary broad¬ 
band PSD, correctly accounting for the intrinsic scatter in the pow¬ 
ers (i.e. equation^, and is computationally very efficient. The ef¬ 
fects of red noise leak can be accounted for by extending the PSD 
model to much lower frequencies than required, generating light 
curves much longer than the observed light curve and using only 
a segment of the required length. Data simulated in such a fashion 
will include power on timescales much longer than those probed by 
each segment alone. 

The fitting procedure is as follows. 


(i) A model PSD is chosen and 500 random light curves are gen- 
eratec^ using the algorithm of Timmer & Kbnig (1995). The simu¬ 
lated light curves were « 26 times longer duration than the orbital 
light curves in order that red noise leak is accounted for. 

(ii) From each of the simulated light curves a section from the 
middle is resampled to match the data. Specifically, from each 
(long) simulated light curve three consecutive light curves were ex¬ 
tracted, to match the three consecutive revolutions, using exactly 
the same sampling as the real data. The logarithmically binned pe¬ 
riodogram is computed in exactly the same way as for the real data 
(Section 4.3). 

(iii) The 500 binned periodograms are then averaged to produce 
the ADM - this represents the periodogram after being distorted by 
the light curve sampling (folded through the window function). 

(iv) The constant Poisson noise level is added and the compari¬ 
son between the ADM and the data is then made by calculating the 

of the fit. 

(v) The normalisation of the ADM is then adjusted to minimise 
the (keeping the Poisson noise level fixed). 

(vi) The x^ of the fit is then evaluated over a grid of PSD models 
and the minimum, corresponding to the best-fitting model, is found. 


This method differs slightly from the method of Uttley et al. 
(2002). Firstly, because the XMM-Newlon light curves are contin¬ 
uous the effects of aliasing are negligible, thus the simulated light 
curves were generated with the same time resolution as the data 
(100 s). The sparsely sampled RXTE data used by Uttley et al. 
(2002) meant that aliasing was a significant factor and needed to 
be accounted for in the simulations. The second major difference 
is that Uttley et al. (2002) calculated errors on the ADM (based 
on the scatter in the simulated periodograms), and these were used 
to estimate the goodness of fit. In the present analysis there are 
enough data points to apply the Papadakis & Lawrence (1993) bin¬ 
ning (with N = 20) and estimate the error on the observed peri¬ 
odogram directly from the data, for use with x^-fitting. 

The x^ fitting employed in this method relies on the data be¬ 
ing binned sufficiently. It would in principle be possible to make 
use of the maximum frequency resolution offered by the unbinned 
periodogram by choosing a different goodness-of-fit estimator (e.g. 
deriving the rejection probability directly from the Monte Carlo tri¬ 
als, as in Uttley et al. 2002). However, the results (the best-fitting 
parameters) are model-dependent in any case and so using a differ¬ 
ent method is unlikely to yield any further information. 


4.3 The periodogram of MCG-6-30-15 

The periodogram of the full-hand data was computed by calculat¬ 
ing the periodogram for each revolution of data (following equa¬ 
tion The light curves from each revolution, when taken sepa¬ 
rately, are perfectly evenly sampled and uninterrupted. This means 
that for each light curve the periodogram estimates are independent 
at each Fourier frequency, and they are also independent between 
the periodograms from the three light curves. The periodograms 
for the three light curves were therefore combined to produce a 
single, binned periodogram by sorting the periodogram points by 
frequency and then applying the logarithmic binning (Papadakis & 
Lawrence 1993). 

® The number of simulated light curves was chosen as a compromise be¬ 
tween the need for an accurate ADM, which requires a large number of 
simulations, and processing time. With 500 simulations, the mean Itr scat¬ 
ter in each ADM frequency bin is 1 per cent. 
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Figure 4. High frequency periodograms estimated using 10 s resolution, 
full-band light curves. At high frequencies the periodograms are dominated 
by the flat (white) Poisson noise spectrum. The solid lines show the ex¬ 
pected Poisson noise levels. 



10 10 ® 10 10 ^ 
Frequency (Hz) 


Figure 5. Periodograms for each energy band calculated using three revo¬ 
lutions of data. The solid lines show the expected Poisson noise levels. 


Furthermore, no additional windowing was applied to the data 
(see Section 13.4 of Press et al. 1992 and section 7.5 of Priest¬ 
ley 1981). Windowing is often used to suppress the leak of power 
from one frequency to another, as it suppresses the side-lobes of 
the spectral window function W{f) and therefore the amount of 
power displaced to distant frequency bins. However, in exchange 
for this the central lobe of the spectral window function is broad¬ 
ened, meaning more power is lost to neighbouring frequency bins. 
This transfer of power between consecutive frequency bins means 
they are no longer independently distributed, which is required by 
the binning method applied above. 

In order to accurately estimate the Poisson noise level, the pe¬ 
riodograms were first calculated using the full-band light curves 
with 10 s resolution. This allows the periodogram to extend up to 
high frequencies, at which the variability is dominated by Poisson 
noise. Fig. Q shows the periodograms from each of the three orbits 
(binned by = 50) and the combined periodogram (binned by 
N — 100 to improve the signal-to-noise). The red noise PSD of the 
source can be seen below ~ 3 x 10“® Hz, above which the (white) 
Poisson noise spectrum dominates. Also marked are the predicted 
noise levels assuming just Poisson noise in the data. These show 
that the predicted noise level is clearly an accurate estimate of the 
true noise level in the data. In the following analysis only the 100 s 
resolution light curves are used, since above the Nyquist frequency 
(5 X 10“® Hz) the periodogram is dominated by the Poisson noise 
and not the source variability, and the Poisson noise level is as¬ 
sumed to be at the expected value. 

Fig. ^ shows the raw periodograms (binned hy N = 20), cal¬ 
culated by combining the periodograms from each revolution. The 
red noise PSD of the source is clearly detected above the Poisson 
noise level at frequencies below ~few 10“® Hz in all bands. In the 
following sections these “dirty” periodograms were fitted using the 
Monte Carlo procedure in order to derive the properties of the PSD 
in a robust manner. 

4.4 Results 

4.4.1 The full-band PSD of MCG-6-30-15 

As expected, the periodograms shown in Fig. show an approx¬ 
imately power-law spectrum for the intrinsic (source) variability. 


Table 3. Results of fitting the full-band periodogram with various trial mod¬ 
els. The models are discussed in the text and shown in Fig. Model 1 is a 
simple power-law, model 2 is a power-law breaking to a slope of ai^vr = 1 
at low frequencies, model 3 is a power-law with a break to aio^ = 0. model 
4 is a power-law with a smooth transition to a slope of = 1 at low 
frequencies and model 5 is an exponentially cut-off power-law. Errors on 
the model parameters correspond to a 90 per cent confidence level for one 
interesting parameter (i.e. a Ax^ = 2.7 criterion). 


Model 

Slope a 

Break freq. /br (10 ® Hz) 

X^/dof 

1 

2.10 

— 

160.8/79 

2 

2.50 ±0.15 

lotF 

87.2/78 

3 

2.45 ±0.15 


89.0/78 

4 

2.65 ±0.20 

13^6 

84.1/78 

5 

1.8±0.I 

100 ± 30 

83.4/78 


Various trial models were fitted to the data (the models are illus¬ 
trated in Fig. ^), using the procedure outlined above to constrain 
the PSD. The first trial model was a power-law (model 1): 

v{f) = Nr‘^. ( 11 ) 

This model has two free parameters, namely the slope (a) and nor¬ 
malisation (N) of the power-law (the additional constant Poisson 
noise level, Pjv, is kept fixed). As summarised in Table ^ this 
model provided an unaccertable fit to the data (> 99.99 per cent 
rejection probability). Fig. M shows the variation in with a (for 
79 degrees of freedom, dof). The of ths fit changes little for 
PSD slopes steeper than 2. This is an effect of the red noise leak. 
If the PSD continues to very low frequencies as a steep power-law 
(a > 2), a significant amount of power will leak into the observed 
periodogram in the form of an a = 2 power-law. This red noise 
leak contribution will cause the observed periodogram to have a 
slope of a « 2 even when the true underlying PSD slope is much 
steeper. Thus, if the power-law PSD continues unbroken to low fre¬ 
quencies, it is difficult to distinguish slopes steeper than 2 in the 
fitting. 

Fig. 1^ shows the best fitting power-law model, with a slope of 
a = 2.1. This steep slope implies that the PSD must break at low 
frequencies or the integrated power would diverge. Therefore PSD 
models including a flattening at low frequencies were explored in 
an attempt to constrain the frequency of the break. The residuals 
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Figure 6. Trial PSD models used to fit the full-band data. 
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Figure 9. The full-band data and (folded) broken power-law model (model 
2 ). 



Slope a 

Figure 7. Fit statistic against power-law index for a simple power-law PSD 
model (model 1, with 79 dof). 



Frequency (Hz) 

Figure 8. The full-band data (points with errors) and best-fitting single 
power-law model (model 1; solid line). The residuals from the fit are shown 
in the lower panel. 


from the power-law fit show that the model is too flat above 2 x 
10~ Hz and too steep at lower frequencies, suggesting the PSD 
may break in the observed frequency range. 

Indeed, based on long timescale RXTE monitoring, Uttley et 
al. (2002) found the best fitting PSD model for MCG-6-30-15, in 
the 2-10 keV band, was a broken power-law. When the slope below 
the break was fixed at aiow = 1 the slope above the break was 
found to bea = 2.0±0.3, with abreakfrequency in therange 1.3x 
10~® to 1.0 X 10~'* Hz (90 per cent confidence limits). Therefore 
the next trial model fitted was a broken power-law (model 2) 



/ ^ /br 
/ > /br, 


( 12 ) 


with a break from a slope of aiow = 1 to a. The free param¬ 
eters were the high frequency slope (a), the normalisation (N), 
and the frequency of the break (/br)- The fit parameters are given 
in Table The model provided an acceptable fit to the data (77 
per cent rejection probability) and the improvement in the fit was 
Ax^ ~ 73.6 for the addition of one free parameter (the break fre¬ 
quency), which, according to the F-test, is significant at > 99.99 
per cent confidence. Fig. shows the confidence contours for the 
model pMameters and the fit residuals using this model are shown 
in Fig. U 

Although there are only two periodogram points at frequen¬ 
cies below the break, the detection of the break is nevertheless 
highly significant. This is because the high frequency PSD slope is 
observed to be steeper than 2, implying that the effect of red noise 
leak is not strong. The steep PSD observed at high frequencies can¬ 
not continue below the lowest observed frequencies; if the PSD 
remained steep to low frequencies the contribution from red noise 
leak would cause the periodogram to resemble a a = 2 power- 
law over the whole observed frequency range. The steep PSD at 
high frequencies must therefore break to a flatter slope close to the 
lowest frequencies probed in order that the high frequency peri¬ 
odogram have a steep slope. 

The slope of the PSD below the break is not well-constrained 
from the XMM-Newlon data however. Fitting the data with a PSD 
model that breaks from a slope of aiow = 0 (model 3) gives a 
similarW acceptable fit but with a slightly lower break frequency 
(Table 0, Fig. MTl). As these two models both give good fits it is not 
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Slope a 

Figure 10. Contour plot showing the parameters of model 2. The contours 
represent 68.3, 90 and 99 per cent confidence levels. 


o 
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Slope a. 

Figure 11. As Fig. [l^ except for model 3. 


possible to distinguish between them based on these data alone. 
However, it is interesting to note that the break frequency found 
from the RXTE data is consistent with the break frequency found 
from the XMM-Newton data when the PSD is assumed to break to 
Qiow = 1- The break frequencies found from fitting the model as¬ 
suming a break to aiow = 0 are inconsistent between RXTE (/br = 
5.112:6 X 10“® Hz) wd XMM-Newton (/br = Slj x 10"® Hz). 

The exact shape of the break is not well defined from these 
data. As an alternative to a sharply broken power-law, a model 
representing a power-law with a smooth break from a slope of 
Qiow = 1 to a was tested (model 4): 


VU) = Nf-^ 



-(oc-l)/2 


(13) 


This has three freytarameters, N, a and /br, and provided an ac¬ 
ceptable fit (Table Fig. [I^, comparable to the models including a 
sharp break, and the break frequency was consistent with that found 
for model 2. 

The final model tested was an exponentially cut-off power-law 
(model 5): 



Slope a 

Figure 12. As Fig. ^except for model 4 



Slope a 

Figure 13. As Fig. [I^except for model 5 


V{f) = (14) 

Again this model has three free parameters, N, a and /br, and pro¬ 
vided a reasonable fit to the data (Table Fig. |^. The best-fitting 
slope for this model is much flatter than for the previous models, 
and the cut-off frequency is much higher than the break frequencies 
in the broken power-law models. That this model provides a good 
fit to the data is perhaps suggesting that the PSD of MCG-6-30-15 
steepens further at high frequencies. 


4.4.2 Energy dependence of the PSD 

On the basis of the XMM-Newton data alone, PSD models 2 through 
5 all provide reasonable fits to the full-band data. In the following 
analysis, model 2 is fitted to the data from the soft, medium and 
hard energy bands, in order to provide a uniform parameterisation 
of the PSD at different energies. The residuals are shown in Fig. 
and the results from these fits are shown in Fig. and Table ^ 
The fits to the soft and hard band data using model 2 are some¬ 
what worse than the fit to the full band data. However, as can been 
seen from Fig. [^, there are no obvious, strong, systematic resid¬ 
uals. Refitting these data with model 4 provided very similar fits 
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Figure 14. Residuals from fits to the soft, medium and hard band data using 
the broken power-law model (model 2). 



Figure 15. Contour plot showing the parameters of the broken power-law 
model (assuming a slope of Qiow = 1 below the break) for the three en¬ 
ergy sub-bands. The contours represent 68.3, 90 and 99 per cent confidence 
levels. 


and refitting using model 5 gave a slightly better fit to the hard 
band data but slightly worse fits to the soft and medium bands. The 
possibility remains that there is some marginally significant struc¬ 
ture in the PSDs at high frequencies (around ~ 4 x 10“® Hz in 
Fig. [^, not accounted for in these simple continuum models, as 
was also suggested for the high frequency PSD of NGC 7469 by 
Nandra & Papadakis (2001). These points notwithstanding, model 
2 was used to parameterise, in a simple and uniform fashion, the 
PSD as a function of energy. 

The break frequencies as determined in each energy band are 
consistent with one another. However, the slope of the PSD above 
the break shows significant changes between the energy bands, with 
the hardest band showing the flattest slope. Table ^ also gives the 
total variability powers, in TVar form, derived from each band (cf. 
Table [^. For each band the total power was estimated by integrat¬ 
ing binned periodogram (after subtracting off the constant Poisson 
noise level) to give TVar,i, and also by integrating the PSD model 
(excluding the Poisson noise component) over the frequency range 
10“® to 5 X 10“® Hz to give T’var,2. As was noted in Section]^ the 
medium band contains the highest total variability power. 


Table 4. Results of fitting the periodograms derived from the four energy 
bands with the broken power-law model (breaking from a slope of a to 
slope of = 1 at low frequencies). Errors on the model parameters 

correspond to a 90 per cent confidence level for one interesting parameter 
(i.e. a Ax^ = 2.7 criterion). Also given are the total powers calculated 
by integrating the binned periodogram (TVar,!) and by integrating the best¬ 
fitting PSD model (T’var,2)- 


Band 

Slope 

a 

X 

1 

o 

xVdof 

-^var,l 

(%) 

-^var,2 

(%) 

full 

2.50 ±0.15 

10+1° 

87.2/78 

26.0 

26.3 

soft 

2.55 ±0.20 

13t^ 

108.4/78 

23.8 

23.8 

medium 

2.30 ±0.15 


89.2/78 

29.3 

29.1 

hal'd 

2.10 ±0.15 

10^1° 

103.5/78 

25.4 

24.0 


4.4.3 Low frequency PSD 

The binned periodograms used above reached down to 3 x 10“® Hz 
in the lowest frequency bin. It would in principle be possible to treat 
the three consecutive light curves as a single, longer light curve 
(top panel of Fig. |i|) and use this to reach even lower frequencies 
in the power spectrum. However, this three revolution light curve 
is no longer uninterrupted, due to the gaps between rev olutio ns. 
This means that the standard DFT (as discussed in Section |4.1.l| ) is 
no longer suitable. There do exist periodogram estimators designed 
to deal with data containing gaps, such as the Lomb-Scargle pe¬ 
riodogram (Lomb 1976; Scargle 1982) but there are drawbacks to 
using such techniques. Firstly, the gaps in the sampling mean the 
periodogram estimates at each frequency are no longer independent 
of one another (Scargle 1982), so the binning and error estimation 
used above is not strictly valid. Secondly, the more complicated 
window function introduces much greater distortion on the low fre¬ 
quency periodogram, which would severely limit the amount of in¬ 
formation that could be recovered from the lowest frequencies. For 
these reasons the results presented in this paper are based entirely 
on the more robust, standard DFT method, but for completeness the 
Lomb-Scargle periodogram (computed using the algorithm of Press 
& Rybicki 1989) of the full-band light curve is shown in Fig. |I| 
Down to low frequencies this appears broadly consistent with the 
results presented above. 


4.4.4 Stationarity 

The fundamental assumption underpinning the above analysis is 
that the variability process is stationary. A stationary process is one 
whose statistical properties do not depend on time (see section 1.3 
of Bendat & Piersol 1986 for more on stationarity in random data). 
Specifically, the assumption made here is that the PSD is station¬ 
ary, i.e. that the light curves are all realisations of the same process. 
The rms-flux correlation found in Section shows that the vari¬ 
ability of MCG-6-30-15 is in some sense non-stationary (the vari¬ 
ance does change), but in a “well-behaved” fashion. Indeed, the 
linear correlation between arms and x means that the average TVar 
is constant. In other words, the variations, when normalised by the 
local flux, appear stationary. 

The method outlined in Appendix A of Papadakis & Lawrence 
(1995) was used as a simple check of the validity of this assump¬ 
tion. If the PSDV(f) is constant then the (unbinned) periodograms 
obtained from each revolution should be identical except for the 
scatter dictated by equation]^ Papadakis & Lawrence (1995) define 
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Figure 16. Unbinned Lomb-Scargle periodogram of the full-ban d ligh t 
curve. The large scatter in the periodogram is expected (see Section hi .3 ). 
The dashed line illustrates the broken power-law model used in Section [4^ . 
Note that the model has not been fitted to the data and that the Lomb-Scargle 
periodo gram is normalised differently from the periodogram discussed in 
Section pl.l.l[ 


the statistic S from the logarithm of the ratio of two periodograms, 
summed over a range of frequencies. 

C. 1 log[A(/i)] -log[-P2(/0] 

VN ^ ^va.T:{log[Pi{fi)]} + var{log[P 2 (/i)]} ’ 

where Pi{fi) and P 2 {fi) are the two periodograms, and 
var{log[Pi, 2 (/i)]} = 0.31 (Papadakis & Lawrence 1993). 

If the underlying PSD is the same in the two periodograms 
then S should be normally distributed with a mean of zero and a 
variance of unity. 

A normalised periodogram was computed for the first 80 ksec 
of each of the three revolutions (the light curves were clipped to 
equal duration so that their periodograms had identical Fourier fre¬ 
quencies) and S was measured for frequencies below 10“^ Hz, 
where the source variability dominates over the Poisson noise. 
Comparing rev. 301 and 302 gave S = 1.18, comparing rev. 301 
and 303 gave S = —0.26 and comparing rev. 302 and 303 gave 
S = —1.44. These three estimates are all within 2a of the expected 
value assuming stationarity in the data. Therefore, to first-order at 
least, the assumption of stationarity in these data seems a reason¬ 
able one. 


5 CROSS SPECTRAL PROPERTIES 
5.1 The cross spectrum 

The PSD is often called the auto spectrum because it is the Fourier 
transform of the light curve, multiplied by itself (equation ^). The 
cross spectrum is a related tool used for comparing the properties of 
two simultaneous light curves, such as from different energy bands, 
say x{t) and y{t). The product of the Fourier transforms of the two 
light curves gives the cross spectrum (compare with equation 

C{f)=XVf)Y{f). (16) 

In the same way that the PSD (auto spectrum) and the ACF are 
Fourier pairs, the cross spectrum and cross correlation function 


(CCF) are also Fourier pairs. Thus, mathematically speaking, the 
cross spectrum contains the same information as the CCF. 

The cross spectrum is complex valued and can be represented 
in the complex plane by an amplitude and a phase. The com¬ 
plex Fourier transforms of the two light curves can be written as 
X{f) = and Y{f) = and the cross 

spectrum can be written as: 

c{f) = (17) 


The phase component of the cross spectrum represents the 
phase difference between the two light curves (which corresponds 
to a time delay at that frequency) and is discussed more fully in 
Section 5.3. The squared magnitude of the cross spectrum is used 


to define the “coherence” of the two light curves. 


5.2 Coherence 

5.2.1 The meaning of coherence 

The coherence, 7^(/), is a real-valued function and is essentially 
the squared cross spectrum normalised by the auto spectra (PSDs) 
of the two light curves. 


(|v(/)P)(|y(/)P)- 


(18) 


Here the angled brackets represent an averaging over an ensemble 
of realisations of the processes. Coherence differs from the peri¬ 
odogram and time lags in the sense that it is only meaningful to talk 
of the coherence from an ensemble of independent measurements 
(either from se par ate light curve segments or consecutive frequen¬ 
cies). Equation!^ defines the linear correlation coefficient of X{f) 
and Y (/) (see Section 9.1 of Priestley 1981)^ The coherence can 
thus be thought of as a measure of the degree of linear correlation 
between the two light curves as a function of Fourier frequency, 
and is related to the amplitude of the CCF. 

In order to illustrate coherence, consider two light curves that 
are related by a transfer function t/>(r), e.g. 


yf) = 


ljj{t — T)x{T)dT = ® x{f). 


(19) 


Then from the convolution theorem 


Y{f) = 'I>(/)A(/), (20) 

where 4'(/) is the Fourier transform of V’(l)- Therefore: 

= Y^if)Y{f) = |v[/(/)|2|X(/)|^ (21) 

and 

C{f) = X*{f)Y{f) = X*{m{f)X{f) = ^{f)\X{ff. ( 22 ) 
The squared magnitude of the cross spectrum then becomes: 
\C{f)\^ = \X{ft\^(f)\\ (23) 

And so 

\c{fr ^ \x{m^{.f)\^ ^ 

\x{mY{fr \x{m<f{mxif)\^ • 


^ The function 7 ^(/) defined by equation [l^ is often referred to as the 
squared coherence (Priestley 1981; Bendat & Piersol 1986) and it gives the 
square of the linear correlation coefficient. This function is simply refeiTed 
to as the coherence in this paper. 
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Therefore the coherence is 7 ^(/) = 1. 

This means that if the two light curves are related by a single, 
linear transfer function (e.g. a delay or a smoothing) then they will 
have unity coherence. In other words, if knowledge of one light 
curve can in principle be used to predict the other, the two are said 
to be perfectly coherent. If, on the other hand, the light curves con¬ 
tain contributions from completely independent processes then the 
coherence will be less than unity. 

Another way of thinking about coherence is in terms of equa¬ 
tion [^. If the light curves are related by some simple transform 
then, at a frequency /, there will be a fixed delay between the light 
curves, i.e. the lag will be the same in each realisation of the pro¬ 
cess. If the phase difference is constant then the value of the numer¬ 
ator of equation (the squared mean of equation is equal to 
the denominator and hence the the coherence is unity. The interpre¬ 
tation of coherence is illustrated more fully in Vaughan & Nowak 
(1997) and Nowak et al. (1999a). 


5.2.2 Estimating the coherence 


In an analogous manner to the PSD, the cross spectrum is estimated 
using the cross periodogram. If the two discretely sampled light 
curves, x(ti) and y(ti), have DFTs and Y(fj) then the cross 

periodogram is given by the discrete form of equation]^, namely: 

C{fj) = X*ifj)Y{fj) = Rcifj) + ilcifj), (25) 

where Rcifj) and Icifj) are the real and imaginary components 
of the complex cross spectrum: 

Rcifj) = Rxifj)RYifj) + Ixifj)lYifj) 

Icifj) = Rxifj)lYifj) - Ixifj)RYifj), (26) 


and Rxifj) and Ixifj) are the real and imaginary components 
of the DFT of the light curve xitj) as defined by equation p . The 
coherence is estimated using the discrete form of equation MS : 


iHfj) = 


(Rcifj))^ + jlcifj))^ 

i\Xifj)?)i\Yifj)\^) 


(27) 


The angled brackets here represent an averaging over light curve 
segments and/or consecutive frequencies. In practice however, the 
coherence will be compromised by photon noise, which is indepen¬ 
dent between the two light curves and so will lead to an apparently 
reduced coherence. As the Poisson noise power is known, this ef¬ 
fect can be removed using the recipe outlined in Vaughan & Nowak 
(1997). In the following analysis the noise-corrected coherence 
and its uncertainty were estimated using equation 8 of Vaughan & 
Nowak (1997). The coherence was calculated only for frequencies 
below 10~® Hz, above this frequency the Poisson noise power be¬ 
comes comparable to the intrinsic source variability and equation 8 
of Vaughan & Nowak (1997) becomes increasingly inaccurate. 

As mentioned above, the coherence needs to be averaged in 
some fashion. For the present analysis the cross and auto spectra 
were calculated separately for each revolution. The data from the 
three revolutions were then combined by sorting in order of fre¬ 
quency and then averaged in bins of = 40 frequency points. 


5.2.3 Coherence properties of MCG-6-30-15 

The coherence function was calculated for the three possible com¬ 
binations of the three energy bands (soft/medium, soft/hard and 
medium/hard). This was repeated using data from MOSl and 
MOS2 as well as the pn and the results were consistent between 



Frequency (Hz) 

Figure 17. Coherence functions between light curves in the three energy 
bands. The different symbols mark the results derived using light curves 
from the three EPIC detectors (squares mark pn, circles mark MOS1 and 
stars mark MOS2 data). The dotted line shows the function 7 ^(/) = 
exp(—//8 X 10“'* Hz) discussed in the text. 
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Figure 18. Coherence functions between light curves from different detec¬ 
tors in the same energy bands (squares mai'k pn vs. MOSl, circles mark pn 
vs. MOS2 and stars mark MOSl vs. MOS2). As these were calculated be¬ 
tween light curves from different detectors but identical energy ranges the 
intrinsic coherence is unity and deviations from unity are a result of Poisson 
noise. The dotted line shows the function plotted in Fig. ^ representing the 
loss of coherence observed between the energy bands. Clearly this function 
shows a much more pronounced loss of coherence than that resulting from 
Poisson noise. 


the detectors. (The coherence was also estimated individually for 
each revolution of data and the results were consistent between the 
three revolutions.) The resulting coherence functions are shown in 
Fig. 0 At low frequencies the coherence is close to unity for all 
energy bands. At higher frequencies the coherence becomes much 
lower, particularly between the soft and hard bands. For compari¬ 
son, Fig. § shows the function 7 ^(/) = exp(—//8 x 10 * Hz) 
compared to the coherence between soft and hard bands. 

This apparent loss of coherence at high frequencies could be 
an effect of Poisson noise in the data, since at high frequencies the 
intrinsic source power is weaker compared to the noise power. In 
order to test this possibility the coherence was evaluated by com¬ 
paring light curves in identical energy ranges from the three dif¬ 
ferent detectors. The light curves extracted from the three detec- 
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tors should be identical, in a given energy band, except for Poisson 
noise. Therefore the coherence should be unity and any deviation 
from unity will be a result of Poisson noise. As can been seen from 
Fig. 0 the coherence is close to unity until the highest frequen¬ 
cies, and the deviation from unity coherence at high frequencies is 
much smaller than the observed loss of coherence between the dif¬ 
ferent energy bands. This, and the fact that the loss of coherence is 
repeated in all three detectors suggests that Poisson noise is not se¬ 
riously biasing these results, and that the error bars are reasonable 
estimates of the true uncertainties. 


5.2.4 Simulations 


Simulated data were used as an additional test of the accuracy of the 
coherence estimation. Two simulated light curves were constructed, 
representing light curves in two different bands, and the coherence 
was calculated for them exactly as for the real data. These simu¬ 
lations were used as a further test of the effect of Poisson noise 
and to measure the effect of different PSD slopes on the coherence 
function. 

The effect of Poisson noise was tested by simulating coherent 
data in two bands using the same PSD. The PSD model used was 
taken from the fit to the soft band data using model 2 (Table ^). Two 
long light curves were generated using this PSD and an identical 
random number sequence. The two light curves, having identical 
PSDs as well as Fourier phases and amplitudes, were thus identical 
and therefore have unity coherence at all frequencies. 

The simulated light curves for the two bands were resampled 
to match the window function of the real data, i.e. three consecutive 
light curves were extracted from each to band match the three rev¬ 
olutions of data. The light curves for the first band were rescaled 
to match the average count rate of the soft band light curve and 
the light curves from the second band were rescaled to match the 
count rate of the hard band data. (The simulated light curves were 
produced using the absolute PSD normalisations found from fitting 
the data. Thus the absolute variability amplitude of the rescaled 
simulated data is as expected based on the real data.) Poisson noise 
was then added independently to the simulated light curves (i.e. the 
number of counts in each time bin was randomised according to 
the Poisson distribution). These two sets of simulated data repre¬ 
sent soft and hard band light curves in terms of their Poisson noise, 
but intrinsically are perfectly coherent and have the same PSD. 

One hundred sets of simulated data were produced and for 
each set the coherence and its error were calculated as discussed in 


Section 5.2.2 (i.e. using the recipe of Vaughan & Nowak 1997). The 


average and standard deviation of the coherence estimates from the 
100 sets of simulated data are shown in the top panel of Fig. 
Also shown are the average error bars. The coherence is clearly 
very close to unity over the frequency range examined, and the 
standard deviation of the estimates is close to the mean error bar, 
indicating that the Vaughan & Nowak (1997) method did accurately 
recover the intrinsic coherence from the effect of the Poisson noise. 
(Note that these error bars are slightly smaller than those derived 
from the real data because the uncertainty is itself a function of the 
intrinsic coherence, which is lower in the real data.) 

The tests above used intrinsically c oheren t data with identi¬ 
cal PSD shapes. The PSDs fits of Section 4.4.2 show that the PSD 


shape is energy dependent. In order to confirm that the coherence 
is not biased because of the PSD shape being different in different 
energy bands, the experiment was repeated using PSDs appropriate 
for the two bands. Specifically, the best-fitting PSD models for the 
soft and hard bands (listed in Table ^ were used to generate simu- 



Frequency (Hz) 

Figure 19. Average coherence function measured from (perfectly coherent) 
simulated data. The top panel shows results using light curves with identical 
PSDs, the bottom panel shows results when different PSDs are assumed for 
the two light curves. The error bars on the data points represent the stan¬ 
dard deviation of the coherence estimates, and the lines without data points 
represent the average of the estimated error bars. Clearly the coherence is 
close to unity and the error estimates are close to the expected scatter. 


lated soft and hard bands light curves, but again the two sets of data 
used an identical random number sequence. Thus the simulated soft 
and hard band light curves again have unity intrinsic coherence but 
different PSDs. As before, 100 sets of simulated data were used 
and the average coherence and error were calculated. The results 
are shown in the bottom panel of Fig. The average coherence is 
again very close to unity, indicating the different PSD slopes in the 
different bands did not bias the results. 

The loss of coherence observed in MCG-6-30-15 is thus not 
an artifact of Poisson noise in the data, nor the change in PSD slope 
with energy, and is repeated in each of the three EPIC cameras. 
As such, in the discussion below, it is assumed that this loss of 
coherence is intrinsic to the continuum variations of the source. 


5.3 Time lags 


5.3.1 Estimating the time lags 


As mentioned in Section 5.1, the cross spectrum also contains in¬ 
formation about the phase differences between the two light curves 
(equation |^. This information is similar to the delays measured 
by the CCF, except measures the delay as a function of Fourier fre¬ 
quency. The phase lag is measured from the phase of the averaged 
cross periodogram: 


=arg{(C(/j))} =tan ^ | {Rc{f\ } ’ 

where C{fj) is calculated using equation ^ and the angled brack¬ 
ets denote an averaging over consecutive frequency points and/or 
light curves. This phase difference can then be converted into a 
time delay at Fourier frequency /: 


T{fj) = <lt{fj)/27tfj. (29) 

Only if the coherence is high at a given frequency is the time lag 
meaningful. If the two light curves are not coherent then the indi¬ 
vidual phase estimates from the two light curves are poorly corre¬ 
lated and so the phase difference of the average cross spectrum has 
no useful interpretation. 
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Figure 20. Estimate of the sensitivity level of time lag measurements. The 
dotted line shows the uncertainty in the time lag due to Poisson noise and the 
solid line represents the appr oxim ate time lags measured, as described by 
the function given in Section 5.3.3 For frequencies above ^ 6 X 10“"^ Hz 
the time lags will not be recoverable from the Poisson noise. 


Figure 21. Time lags between light curves in the three energy bands. The 
symbols are as in Fig. |^. A positive lag indicates the first band leads the 
second. The dotted line shows the function t(/) = 0.013/“^ discussed in 
the text. 


For the present analysis, the time lag was calculated using 
equations ^and^, with the averaging of the cross spectrum per¬ 
formed in bins of Af = 40 frequency points. The error on the 
time lag estimate was computed using equation 16 of Nowak et al. 
(1999a). See section 9.1.3 of Bendat & Piersol (1986) for a deriva¬ 
tion of this equation. 


5.3.2 Poisson noise ejfects 

The effect of Poisson noise is to add an independent, random com¬ 
ponent to the phases of the two light curves and hence randomise 
the phase difference. This means that small time lags are difficult 
to detect due to photon noise. The recipe outlined in Nowak et al. 
(1999a) (see their Section 4.2) was used to estimate the effective 
noise limit, below which it is difficult to detect reliable lags. The 
estimated error on the phase lag is computed using their equation 
14: 

This uncertainty was calculated using values appropriate for the 
hard band (PSD normalisation, slope and noise level Pn from Sec¬ 
tion [l.4t a nd coherence estimated using the function given in Sec¬ 
tion p(2^. This estimated uncertainty, which should be considered 
merely a lower limit on uncertainty introduc ed by Poisson noise, 
was compared to the function used in Sectionto illustrate the 
observed time lags and the results are shown in Fig. This indi¬ 
cates that above ~ 6x10”'* Flz the lags are too small to be detected, 
i.e. they fall below the sensitivity limit defined by the Poisson noise. 
A further point is that the coherence is significantly below unity at 
higher frequencies (^few 10“* Hz), in which case any measured 
time lag is no longer meaningful. Therefore, from these data, sig¬ 
nificant time lags will only be detectable below ~few 10”* Hz. 


5.3.3 Time lags in MCG-6-30-15 

The time lags were computed (for each EPIC detector) for the three 
combinations of energy bands. The results are shown in Fig. 0 


In each test (three combinations of energy bands and three differ¬ 
ent detectors) there was a significant lag at the lowest frequencies, 
with the softer band leading the harder band. At higher frequencies 
the lag becomes consistent with zero. The lag is small, only about 
~ 1 per cent between soft and hard bands (e.g. on a timescale of 
20,000 s the lag is only ~ 200 s). This is illustrated in Fig. 
which shows the function r(/) = 0.013/”* compared to the mea¬ 
sured time lags between the soft and hard bands. 

The time lags were calculated using light curves for the three 
different detectors in identical energy ranges , to test the possible 
influence of Poisson noise (see Section 5.2.3). These are shown in 
Fig. 1^. In these tests the lag was consistent with zero over most of 
the useful frequency range. This, combined with analysis of Sec¬ 
tion 5.3.2 and the fact that the lags were recovered independently 
from each EPIC camera, suggests that the lags observed between 
different energy bands are not an artifact of Poisson noise (which 
in any event would destroy any time lag signal). 


5.3.4 Simulations 

The simulated data used in Section ^.2.4 to assess the accuracy of 
the coherence function were also used to assess the accuracy of the 
time delay spectrum. 

In the first set of simulations the data in the two bands were 
generated using the same PSD and the same random number se¬ 
quence. In addition to having unity coherence, the pairs of sim¬ 
ulated data also have identical Fourier phases and therefore zero 
intrinsic phase lag between them. The top panel of Fig. ^ shows 
the average of the time delay estimates, their standard deviation 
and the average error estimate from the 100 pairs of sim ulated data 
(including Poisson noise as described in Section 5.2.4). This aver¬ 
age spectrum is close to the zero-delay expectation meaning that 
the Poisson noise did not bias the time delay measurements, and 
the average estimated error is close to the standard deviation of the 
estimates. 

In the second set of simulations the pairs of simulated data 
were generated using different PSDs for the two bands. The bot¬ 
tom panel of Fig. ^ shows the average time delay spectrum, which 
again is as expected if the delay estimate is unbiased. It therefore 
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Figure 22. Time lags between light curves from different detectors in the 
same energy bands. Squares mark pn vs. MOSl, circles mark pn vs. MOS2 
and stars mark MOSl vs. MOS2. As the light curves being compared cover 
identical energy ranges they should differ only as a result of Poisson noise 
as so result in zero lag, any lags obseiwed will be a result of Poisson noise. 
The dotted line shows the function plotted in Fig. Q. As can be seen, the 
observed lags between bands (Fig. Elh are much larger than those resulting 
from Poisson noise. 
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Figure 23. Average time delay measured from simulated data (with zero 
intrinsic delay). The top panel shows results using light curves with identi¬ 
cal PSDs, the bottom panel shows results when different PSDs are assumed 
for the two light curves. The error bars on the data points represent the stan¬ 
dard deviation of the time delay estimates and the lines without data points 
represent the average of the estimated error bars. The phase is close to zero 
and the en'or estimates are close to the expected scatter. 


seems reasonable to assume that the time delays observed in MCG- 
6-30-15 at low frequencies are intrinsic to the source. 


6 DISCUSSION 
6.1 Summary of results 

This paper presents an analysis of the X-ray continuum variabil¬ 
ity of MCG-6-30-15. A strong, linear correlation is found between 
the rms variability amplitude and the source flux (Section 3.2). The 
PSD was examined by fitting model power spectra to the data af¬ 
ter allowing for the distorting effects of sampling using a Monte 



Frequency (Hz) 

Figure 24. Power spectral data from the three energy bands shown in 
fP{f) units (dimensionless). The data have be “unfolded” using model 
2 to better illustrate the shape of the underlying power spectra. The PSD 
models are shown with dotted lines, and the error bars show the data/model 
residuals multiplied by the unfolded PSD model. The bottom right panel 
shows the three model PSDs overlayed. 


Carlo procedure (Section However, the excellent sampling of 
the XMM-Newton light curves means that these biases are in any 
case minimised and so there was a minimum of information lost 
from the PSD. The PSD is well represented by a power-law of slope 
aiow = 1 at low frequencies, breaking to a slope of a « 2.5 at a 
frequency /br ~ 10““^ Hz. However, it should be re-iterated that 
these PSD results are model-dependent, and that parameterising the 
PSD using a different model yields slightly different results. 

The PSD is energy dependent, with the variations in the higher 
energy band showing a flatter PSD slope (above the break fre¬ 
quency) than variations in the softer bands. There is no evidence 
for any energy dependence of the break frequency. Fig. ^ illus¬ 
trates the PSD shape and energy dependence. Fitting the PSD as 
a function of energy revealed tentative hints of ene rgy-dependent, 
high-frequency structure in the PSD (Section 4.4.2). This could be 
due to subtle instrument effects or may represent low-amplitude 
structure intrinsic to the source PSD (see also Nandra & Papadakis 
2001). These timescales (^ 10“® Hz) are comparable to the dy¬ 
namical timescale at innermost stable orbit around the black hole. 
However, further progress on the very high frequency variability of 
relatively faint sources such as Seyfert galaxies must await future 
high-throughput missions such as Constellation-X and XEUS. 

The variations observed in different energy sub-bands are 
highly coherent at low frequencies (long timescales) but the co¬ 
herence falls significantly below unity above ~few 10“"^ Hz, close 
to the break in the PSD (Section 5.2). The high coherence between 
soft and hard bands at low frequencies means that the two bands 
are linked by a single, linear transfer function. The direction and 
magnitude of the time lag shows that the hard photons lag the soft 
photons, on average, by ~ 200 s at a frequency / « 6 x 10“® Hz 
(Section 5.3). It is not possible to reliably constrain the frequency 
dependence of the time lags, but the data seem consistent with an 
approximately r oc dependence of the lag (i.e. constant phase 
difference 4>{f)) similar to that seen in the Seyfert galaxy NGC 
7469 (Papadakis et al. 2001) and in GBHCs (Nowak et al. 1999a; 
Pottschmidt et al. 2000). On shorter timescales the coherence falls 
well below unity implying there is no simple transfer function re¬ 
lating soft and hard bands. 
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6.2 Comparison with previous analyses 


Uttley et al. (2002) presented an analysis of long timescale mon¬ 
itoring of MCG-6-30-15 using RXTE. The PSD was found to be 
best described by a broken power-law with slope of aiow = 1 
at low frequencies, breaking to a = 2.0 ± 0.3 at a frequency 
/br = 0.13 — 1.02 X 10“"* Hz. Using the same PSD model, the 
break frequency found from the high-frequency XMM-Newlon data 
is /br = 0.6 — 2.0 X 10~^ Hz, consistent with the RXTE result. 
The slope of the PSD above the break is different, but this is most 
likely an effect of the energy dependence of the PSD. Indeed, us¬ 
ing the 2-10 keV XMM-Newton data (i.e. the same energy range 
used by Uttley et al. 2002) the slope of the high frequency PSD is 
a Ri 2.10 ± 0.15, which is consistent with the RXTE result. 


The break frequency presented here (and in Uttley et al. 2002) 
is model dependent. This is because a model PSD had to be as¬ 
sumed and fitted to the data to derive the appropriate parameters. 
However, the broken power-law model (with a break to aiow = 1 
at low frequencies) provides consistent results when applied to the 
long timescale RXTE monitoring and the short timescale XMM- 
Newton observation. This, and the similarity of the observed PSD 
to those of GBHCs (see Section 5.3) suggests that the chosen model 
(and therefore the derived break frequency) is reasonable. 


Nowak & Chiang (2000) used RXTE and ASCA data to mea¬ 
sure the PSD of MCG-6-30-15. They found the PSD to break from 
a~0toa« lata frequency ~ 10“® Hz and then break again 
to a « 2 in the range 10“"^ — 10“® Hz, The position of the high 
frequency break, as well as the PSD slopes, are roughly consistent 
with those found from the XMM-Newton data. However, as noted 
by Uttley et al. (2002), the method employed by Nowak & Chiang 
(200 0) did not account for distorting effects of sampling (see Sec¬ 
tion 1.1.4) and underestimated the errors at low frequencies. Once 
these effects are included the low-frequency flattening found by is 
no longer significant (section 6.2 of Uttley et al. 2002). Nowak & 
Chiang (2000) also place a limit on any time delay between soft 
(0.5-2.0 keV) and hard (8-15 keV) bands of r < 2000 s, consis¬ 
tent with the time lag results presented here. 


Hayashida et al. (1998) also estimated the PSD of MCG-6- 
30-15, this time using Ginga data. They found a broken power-law 
to be a good fit to the data. However, these authors also did not 
account for the biases in the periodogram, which will be significant 
due to the patchy sampling of the Ginga light curves. The unusually 
flat low-frequency slope (a = —0.5) found by these authors is 
most likely a result of bias in the low-frequency periodogram, and 
therefore the position of their claimed break frequency should be 
considered with some caution. 


The only other AGN for which detailed cross spectral analy¬ 
sis has been performed is NGC 7469 (Nandra & Papadakis 2001; 
Papadakis et al. 2001). In this object the PSD does not show any 
significant breaks, but does show the same energy dependence as 
found in MCG-6-30-15. The present analysis shows that the en¬ 
ergy dependence of the PSD is due to the slope of the PSD be¬ 
coming flatter at higher energies, rather than a change in the break 
frequency. In addition, NGC 7469 shows time lags between energy 
bands in the same direction (soft leading hard) and with the similar 
relative magnitude (~ 1 per cent) as in MCG-6-30-15. Papadakis 
& Lawrence (1995) showed that in NGC 4051 the harder band 
showed a flatter PSD than the softer band. Recently, M'^ Hardy et 
al. (in prep.) used XMM-Newton data of NGC 4051 and also found 
this energy dependence of the PDS. 


6.3 Comparison with Cygnus X-1 

Since EXOSAT first showed the medium timescale PSDs of Seyfert 
galaxies to have a power-law form, their timing properties have 
been compared to those of GBHCs. Now, with long timescale mon¬ 
itoring available from RXTE (Edelson & Nandra 1999; Nandra & 
Papadakis 2001; Uttley et al. 2002) and short timescale observa¬ 
tions possible with XMM-Newton (this paper) it is possible to re¬ 
assess this comparison. 

The best studied GBHC, Cyg X-1, when in its low/hard state 
shows a PSD with an approximately broken power-law form. The 
PSD breaks from a low frequency slope of a ~ 0 to a « 1 and then 
to Q ~ 2 at high frequencies. The position of the high frequency 
break is typically /br « 3 Hz (Belloni & Hasinger 1990; Nowak et 
al. 1999a). The comparable break in the PSD of MCG-6-30-15 oc¬ 
curs at /br ~ 10“^ Hz. This suggests the characteristic timescales 
in these two systems differ by a factor ~ 3 x 10"*. Under the as¬ 
sumption that the timescales scale linearly with mass of the central 
black hole, and assuming Cyg X-1 to contain a 10 Mq black hole 
(Herrero et al. 1995), then suggests the mass of the black hole in 
MCG-6-30-15 is ~ 3 X 10® Mq (in the range 1.5-5.0 x 10® M© 
using the 90 per cent limits on the break frequency in MCG-6-30- 
15). Assuming a bolometric luminosity of 4 x 10"*® erg s“® for 
MCG-6-30-15 (Reynolds et al. 1997; converting to Hq = 75 km 
s“® Mpc“®) then implies that MCG-6-30-15 is radiating at the 
Eddington luminosity. However, since independent mass estimates 
(such as from reverberation mapping) are not available for MCG- 
6-30-15 it remains possible that this mass estimate is inaccurate if 
the relevant timescales do not scale linearly with black hole mass 
(as assumed above). 

If the frequency of the PSD break is related to the orbital fre¬ 
quency of matter close to where the peak emissivity of the accretion 
disc occurs, rpeak, then it will depend on rpeak/fin, where rin is the 
inner radius of the disc. Eor a standard disc, rpeak/fin ~ 2. This 
will also introduce a dependence on rin/rg (where rg = GM/c? is 
the gravitational radius of the black hole) and thus the spin parame¬ 
ter of the black hole a/M. If the black hole in MCG-6-30-15 has a 
higher spin parameter than that in Cyg X-1 then the mass deduced 
from the above scaling could be increased. 

The other temporal characteristics seen in MCG-6-30-15 are 
remarkably similar to those observed in Cyg X-1, implying there is 
indeed a strong connection between these two types of source. The 
PSD (above the break frequency) is flatter in harder energy bands 
in both MCG-6-30-15 and Cyg X-1 (Nowak et al. 1999a; Lin et al. 
2000). The strong rms-flux correlation first seen in X-ray binaries 
(Uttley & M'^Hardy 2001) is exhibited by MCG-6-30-15 (see Sec- 
tion|3.2|). MCG-6-30-15 also displays a loss of coherence between 
hard and soft bands at high temporal frequencies (i.e. above the 
break in the PSD) very similar to that seen in Cyg X-1 (Nowak et 
al. 1999a) and time lags between these bands that appear similar (in 
direction and relative magnitude) to those observed in Cyg X-l. In 
particular, the time lags measured in MCG-6-30-15 between soft 
and hard bands is ~ 200 s in the lowest frequency bin, where the 
coherence is near unity. In Cyg X-1 the shortest time lags observed 
are ~ 2 X 10“® s (Nowak et al. 1999a; Pottschmidt et al. 2000) at 
frequencies just below where the coherence falls significantly be¬ 
low unity. The ratio of these two time lags is ~ 10®, similar to the 
ratio between the PSD break timescales. 

It should also be noted that it is not yet clear whether Seyfert 
galaxies more closely resemble the low/hard or high/soft states 
of GBHCs (see also the discussion in Uttley et al. 2002). In the 
high/soft state the PSD of Cyg X-1 is steep at high frequencies 
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(a > 2) and flattens to a « 1 at fhi ~ 10 Hz, with no further 
flattening evident at lower frequencies (Cui et al. 1997; Churazov 
& Gilfanov & Revnivtsev 2001; Reig, Papadakis & Kylatis 2002). 
This ratio of characteristic timescales in this case is ~ 10® and the 
corresponding mass estimate for MCG-6-30-15 is ~ 10® M©, im¬ 
plying a luminosity ~ 30 per cent of the Eddington limit. This high 
(but nevertheless sub-Eddington) accretion rate might be expected 
if MCG-6-30-15 is indeed in a state similar to the high/soft state 
seen in Cyg X-1. 

It is interesting to note that the best-fitting linear model to the 
rms-flux correlation observed in MCG-6-30-15 is consistent with 
passing through the zero-flux origin (Section 3.2). The simplest ex¬ 
planation of this is if there is no significant constant component to 
the light curve. A strong constant component (~ 25 per cent of 
the source flux) is observed from the rms-flux correlation in Cyg 
X-1 (Uttley & M'^Hardy 2001), but only when the source is in the 
low/hard state (Gleissner et al. 2002). This is perhaps also indicat¬ 
ing that the variability properties of MCG-6-30-15 are more similar 
to the high/soft state than the low/hard state of Cyg X-1. (The time 
delay properties of Cyg X-1 when in its high/soft state are rather 
similar to those seen in its low/hard state, see e.g. Pottschmidt et al. 
2000 ). 


At higher temporal frequencies the coherence falls well be¬ 
low unity, meaning there is no longer a simple transfer function 
between the soft and hard band emission. As a static corona should 
maintain a transfer function (whatever the physical mechanism), 
on timescales longer than the break in the PSD (^ 10"^ s), where 
the coherence is high, the corona may be effectively static (i.e. it 
maintains its average properties on these timescales). On shorter 
timescales the corona may be dynamic and so the steep PSD above 
the break and the loss of coherence may be due to changes in the 
corona itself. 

The timescale of the break in the PSD from f~^ to a steeper 
slope is 10“* s, which corresponds to a light-crossing distance of 
~ 2 X lO'^rg. The orbital timescale at 40rg around a 10® Mq black 
hole is torh ~ 10® s, while the thermal timescale for a standard, 
thin accretion disc is fth ~ 10“^ s (assuming a viscosity parameter 
a = 0.1), which is of the same order as the timescale of the break 
in the PSD. The viscous timescale in such a disc is many orders 
of magnitude larger than these other timescales (see the discussion 
of accretion disc timescales in section 5.8 of Frank, King & Raine 
1985). However, if that part of the accretion flow responsible for 
modulating the X-ray emission is geometrically thick, the break in 
the PSD could correspond to a viscous timescale (which depends 
on the ratio of disc height H to radius i? as ~ {HjRf)- 


6.4 Physical implications 

The X-ray emission mechanism operating in Seyfert galaxies is 
usually thought to be inverse-Compton scattering of soft photons in 
a hot corona (e.g. Sunyaev & Titarchuk 1980; Haardt & Maraschi 
1991). In the simplest models the harder photons are expected to lag 
behind the softer photons due to the larger number of scatterings re¬ 
quired to produce harder photons, and the delay should be of order 
the light-crossing time of the corona. The longest observed time 
delay in MCG-6-30-15 between soft and hard bands is ~ 200 s. 
Assuming a black hole mass of 10® Mq, this lag corresponds to 
a light-crossing distance of ~ 40rg (rg = GM/c?). Thus the di¬ 
rection and magnitude of the observed lag in MCG-6-30-15 are 
consistent with an origin in a Comptonising corona. However, if 
the lags are frequency dependent (as expected by analogy with Cyg 
X-1 and also seen in NGC 7469; Papadakis et al. 2001) the lags 
at lower temporal frequencies would become much longer than 
expected for a compact corona (see discussion in Nowak et al. 
1999b). Indeed, the simplest such models predict the time delay 
between soft and hard photons to be independent of Fourier fre¬ 
quency (Miyamoto et al. 1988), contrary to what is observed. In 
addition, some models of Compton scattering coronae predict the 
high frequency PSD should be steeper for higher energy photons, 
due to the high-frequency fluctuations being washed out by multi¬ 
ple scatterings (Hua & Titarchuk 1996; Nowak & Vaughan 1996), 
again contrary to the observations. 

Alternatively, the time delay between soft and hard bands 
could be due to the spectral evolution of individual X-ray events. If 
the spectrum of the event becomes harder as the event unfolds this 
can lead to an average lag between soft and hard bands. Such mod¬ 
els have been discussed by Bdttcher & Liang (1998) and Poutanen 
& Fabian (1999). 

Another alternative is that the delay between soft and hard 
bands could be due to a propagation effect, i.e. a “trigger” signal 
reaching the soft emission region before the hard emission region. 
In this scenario the magnitude of the lag is related to the propaga¬ 
tion speed and the distance between soft and hard emission regions. 
The following section describes a model based on this kind of sig¬ 
nal propagation. 


6.5 A phenomenological model 

Various models have been considered to explain the variability 
of GBHCs and Seyfert galaxies. These include: various exam¬ 
ples of “shot noise” (e.g. Terrell 1972; Lehto 1989; Merloni & 
Fabian 2001); rotating “hot spots” on the surface of an accretion 
disc (Abramowicz et al. 1991; Abramowicz 1992; Bao & 0stgaard 
1995); occultation by moving clouds (Abrassart & Czerny 2000); 
self-organised criticality in an accretion disc (Mineshige, Ouchi 
& Nishimori 1994); fluctuations in the accretion rate propagating 
through the disc (Lyubarskii 1997; Churazov et al. 2001; Kotov, 
Churazov & Gilfanov 2001) and magnetohydrodynamic instabili¬ 
ties in the inner accretion disc (Hawley & Krolik 2001). 

Most of these models were developed to explain only one as¬ 
pect of the variations (e.g. its red noise nature). Now, however, it is 
possible to compare the predictions of these models with the well- 
determined timing properties of MCG-6-30-15. A viable physical 
model has to be able to reproduce not only the correct PSD shape 
and timescales but also the other temporal characteristics such as 
the rms-flux correlation, the energy dependence of the PSD, the 
time lags and the coherence. It is not clear whether any of the above 
mechanism can explain all the observations, and indeed some mod¬ 
els (e.g. the obscuration model of Abrassart & Czerny 2000) seem 
to offer little hope of explaining properties such as the coherence 
and time lags. 

A simple phenomenological model has been developed by 
Lyubarskii (1997), Churazov et al. (2001), and Kotov et al. (2001) 
that seems to explain (qualitatively, at least) the temporal properties 
of Cyg X-l. As the timing properties of MCG-6-30-15 discussed 
above are very similar to those of Cyg X-1, it is worth discussing 
this model as a possible explanation for the behaviour of Seyfert 
galaxies. 

In its high/soft state, Cyg X-1 shows a f~^ PSD over several 
decades in frequency (Cui et al. 1997; Churazov et al. 2001; Reig 
et al. 2002), down to timescales far longer than would normally be 
associated with the X-ray emitting region itself (which is assumed 
to be compact). The model put forward by Lyubarskii (1997) to ex¬ 
plain this very broad range in timescale has variations in accretion 



18 Vaughan et al. 


rate occurring over a large range of radii from the central black hole 
and propagating in towards the X-ray emitting region (assumed to 
extend over only a few tens of gravitational radii). 

If the viscosity parameter of the accretion flow is varying ran¬ 
domly and with the same amplitude at all radii, variations in mass 
accretion rate at much smaller radii will have an f~^ PSD. The 
accretion rate at the X-ray emitting region (assumed to be concen¬ 
trated at small radii) is thus being modulated by variations pro¬ 
duced over a large range in radius. If the X-ray emission emerg¬ 
ing from this region is proportional to the local accretion rate then 
the X-ray variations will also have an f~^ PSD extending down to 
low frequencies. Each radius produces variations at a characteris¬ 
tic timescale comparable to the timescale for inward propagation, 
hence long-timescale fluctuations are produced at large radii and 
shorter-timescale fluctuations are produced closer in. This works 
because any variations on timescales shorter than the propagation 
timescale are damped as they move through the disc (see Lyubarskii 
1997, and Churazov et al. 2001). This naturally explains both the 

part of the PSD (i.e. below the high frequency break; Chura¬ 
zov et al. 2001) and the rms-flux correlation (Uttley & M'^Hardy 
2001 ). 

Once these variations reach the inner regions, where accre¬ 
tion energy is released as X-ray emission, variations on shorter 
timescales are further suppressed, leading to a steepening of the 
PSD. Thus the high frequency break in the PSD from to a 
steeper slope (as seen in MCG-6-30-15) is due to suppression of 
fluctuations within the X-ray emitting region, and the timescale of 
the break in the PSD corresponds to the characteristic timescale at 
the outer radius of the X-ray emitting region (Churazov et al. 2001). 
The emissivity as a function of radius defines a Green’s function for 
the incoming f~^ fluctuations. In frequency terms, the radial extent 
of the X-ray producing region makes it act as a “low-pass filter” for 
the incident f~^ fluctuations, resulting in a break in the PSD, to a 
steeper slope, at high frequencies. 

The specific model described by Churazov et al. (2001) has a 
relatively stable, optically thick, geometrically thin accretion disc 
sandwiched by an optically thin, geometrically thick corona. The 
disc produces relatively slowly varying thermal emission (which 
would emerge presumably in the optical/ultraviolet for a Seyfert 
galaxy) while it is the variations within the corona that propagate 
inwards to produce the X-ray variability. The timescale of the high 
frequency break thus corresponds to the accretion timescale (~the 
viscous timescale) of the geometrically thick corona at the outer 
edge of the X-ray emitting region. 

Kotov et al. (2001) extended this model by adding the simple 
condition that the spectrum of the X-ray emitting region is a func¬ 
tion of radius, with the outer regions producing a softer spectrum 
than the inner regions. This model can then reproduce the observed 
time lags between soft and hard emission. Accretion rate variations 
propagating inwards reach the softer emitting region (larger radii) 
before the harder emitting region (smaller radii), explaining why 
the soft emission leads the hard. This can explain the energy de¬ 
pendence of the time lags (as delay increases with the separation 
between energy bands) as well as the frequency dependence (since 
different radii correspond to different characteristic frequencies as 
well as different energy spectra). In addition, this model explains 
the energy dependence of the PSD above the high frequency break. 
The harder emission is produced from smaller radii, which can be 
modulated by high frequency fluctuations, hence the harder emis¬ 
sion will have more high frequency power than the softer emission 
(i.e. the hard band PSD will be flatter), as observed (see Fig. 0- 
Also, as the variations on short timescales are occurring within 


the X-ray emitting region, i.e. are a result of the dynamic corona, 
the coherence between bands will fall off at high frequencies, as 
mentioned in Section 6.4 (see also the discussion in Nandra & Pa- 
padakis 2001). 

Such phenomenological models require that variations occur¬ 
ring at large radii can be fed into much smaller radii and so in¬ 
fluence the PSD. The observed iron line profile of MCG-6-30-15 
from this observation (Fabian et al. 2002) indicates that the disc ex¬ 
tends down to less than 2rg, and, from the steep inferred emissivity 
profile, that much of the X-ray power is generated very close to 
the black hole. This presents a challenge for all such models, given 
that the dynamical timescale of this innermost region is 10-100 s, 
whereas the timescales measured (e.g. the PSD break timescale) 
are orders of magnitude longer. If the f~^ PSD of MCG-6-30- 
15 is similar to that in the high/soft state of Cyg X-1 (Reig et al. 
2002) then the variability timescales will extend to many orders of 
magnitude longer still. Indeed, the total power in f~^ spectra di¬ 
verges so slowly that the PSD could in principle continue as f~^ 
until very low frequencies without the total power becoming large. 
(For example, model 2 from Section 4.4 could continue down to 
the unphysically low frequency ^ '' 
power exceeding F^ai = 1!) 


10 Hz without the integrated 


7 CONCLUSIONS 

This paper presents a detailed time series analysis of the contin¬ 
uum variability of MCG-6-30-15. The PSD shows a break from 
f~^ to/“^ ® at a frequency /br ~ 10“"^ Hz. Comparing the break 
timescale with the analogous timescale seen in Cyg X-1, and as¬ 
suming that these scale linearly with black hole mass, gives a black 
hole mass of ~ 10® Mq for MCG-6-30-15 and suggests it is ac¬ 
creting at a significant fraction of the Eddington limit. The slope 
of the PSD above the break is energy dependent, with the harder 
band showing more high frequency power than the softer band. The 
variations in different energy bands also become incoherent above 
the break in the PSD. This leads us to speculate that on timescales 
shorter than the break timescale, the X-ray emitting corona is dy¬ 
namic, and that harder photons are produced at smaller radii than 
softer photons. At the lowest Fourier frequencies the soft and hard 
bands are highly coherent and the soft band is found to lead the 
hard band by ~ 200 s. The direction and magnitude of this lag is 
consistent with simple Comptonisation models. A phenomenologi¬ 
cal model, originally developed to explain the PSD of GBHCs 
such as Cyg X-1, is discussed and found to be able to reproduce 
(qualitatively) the observed temporal characteristics of MCG-6-30- 
15. Although this model is rather simple and phenomenological it 
can explain the observed temporal characteristics of the X-ray con¬ 
tinuum in MCG-6-30-15. 
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